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Abstract 

A fit is made to the data for the proton structure function up to = 10 GeV^, including the 
real 7]? total cross-section. It is economical and simple, and its form is motivated by physical 
principles. It is extrapolated down to very small values of x. Data for the ratio z^VF^/z/VF|' 
are also fitted. 



A FORTRAN program for the Gt to z/VF|' is available by email on request 



May 1993 



tad@ v2.ph.man.ac.uk 
*pvl@amtp. cam.ac.uk 



The structure function h'W2 for inelastic electron or muon scattering is constrained by gauge 
invariance to vanish linearly with at — 0. Indeed, the total cross-section for real-photon 
scattering is 

(^ilP) = vW2 (1) 

evaluated at Q"^ = 0. For this reason, if for no other, the variation of vW2 with at small 
Q"^ cannot be described by perturbative QCD: it is unsafe to use any perturbative evolution 
equation until is at least so large that vW2 has fully recovered from the need to vanish at 

A few years agot"*^] we parametrised the then-available data for cr(7p) and for i'W2 at small 
together in a very simple form. A rather similar analysis has been repeated more recently!^]. 
Our parametrisation compares quite well with more recent data, but we feel that the time 
has now come to improve on it. This is for two reasons. The first is concerned with the search 
for the Lipatov pomeron in the measurements of vW2 at small x. There is some uncertainty 
about the effects of kinematic constraintsI'^lt'^J. and of shadowing on the solution to the 
Lipatov equation, and so it may not be immediately obvious from future data whether the 
Lipatov pomeron is actually present. To help decide this, it will be necessary to know as 
accurately as possible what is expected from the soft pomeron, which is the pomeron that is 
seen in all the data that have been collected so far[6]. Our new parametrisation allows what 
we believe to be a reliable extrapolation to very small x of the effect of soft pomeron exchange. 
The second reason good fits to small-Q^ data are needed is for calculation of QED radiative 
corrections. These can be huge at HERA and so they must be calculated as accurately as 
possible, using more than one parametrisation of existing data as a check. 

Our aim is to present a parametrisation of the data that is as simple as possible, with rather 
few parameters, and whose component parts are motivated by physical considerations. 

We have shown recently!^] that all total cross-sections may be parametrised as a simple sum 
of two Regge powers 

We fixed the two powers in this expression from pp and pp data, and found that the resulting 
expression works well also for irp and Kp. In the case of 7p, the fit (2) is very similar to 
one we made a few years agot"*^], and it is in agreement with the published data points from 
hera[^]. We show this in figure la. (Notice that although we have specified the two powers 
to high accuracy, the data do not determine them to such accuracy; the pp and pp data from 
which we derived them would allow slightly different values, with corresponding changes to 
the coefficients X and Y - the errors in the parameters are strongly correlated.) 
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According to the parton model, the same Regge powers appear!^] as powers of 1/x in the 
small-x behaviour of i'W2- The parton model is drawn in figure 2, where the amplitude T 
is the amplitude for the emission of a parton of momentum k by the proton. From simple 
kinematics [^], the energy variable for this amplitude is 

{p - kf = -x-\k^ + \4) -k^ + {l- x)ml (3) 

and so is large when x is small. Because T is a hadronic amplitude, its high-energy behaviour 
should involve the same Regge powers of {p — k)"^ as appear as powers of s in (2), and these 
reflect themselves as corresponding powers of 1/x in i'W2- Of course the parton model is only 
the first term in a perturbative-QCD expansion of vW2', the infiuence of the perturbative 
corrections on the Regge powers is not understood and is the reason for the great interest 
in small-x physics at HERA. One possibility is that the Regge powers in (2) survive the 
perturbative corrections and appear in the small-a; behaviour, at least for small, and even 
perhaps moderately large, values of Q^. We explore this possibility in this paper. 

We note that Monte Carlo models are commonly in error by taking both k"^ and k^- to be 
zero. While this may be a good approximation for many purposes, it is not good when x is 
small, because, according to (3), it fails to make (p — kY large and so does not correctly take 
account of the nonperturbative Regge behaviour. 

Our first fit, then, is to the small-x data from We use two simple powers of each 

multiplied by a simple function of that vanishes linearly with as — > and goes to 
1 for large Q^: 

. ^2 X 1-0808 . ^2 X 0.5475 

. ^.-cosos (_^L_) (4a) 

with the constraints that 

^^-1.0808 ^ Q_gQ4 56-0-^^^5 = 1.15 (46) 

so as to retrieve the fit of figure la to the real-photon data when 0. This two-parameter 

fit is compared with the NMC data in figure 3, for the choices A = 0.324 and B = 0.098. 
Our fit has used only the data up to = 10. 

Strictly speaking, because x = Q'^/2i', the ^ limit of (4) is not a sum of powers of s as 
in (2), but rather a sum of powers of 2i/ = s — m%. At the lowest s- value of the data in figure 
la the diff'erence is negligible, though if we continue to smaller s it becomes important. We 
show the curve obtained for a{'yp) from the — limit of (4) in figure lb. Even though 
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we have not used the data below ./s = 6 to determine the fit, it works satisfactorily down 
into the region where the resonances begin to be important. 

If we extrapolate (4) to very small x, we obtain a predicted value just less than 0.8 for ^^^2 
at a; = 10~^, = 10. This prediction remains rather stable under the refinements to our fit 
that we report below. We make refinements to the fit partly in order to extend it to larger 
values of x, as is necessary if it is to be useful for making radiative corrections. But also 
we want to incorporate the contribution from heavy flavours into the analysis, since part of 
the rapid rise with seen in the data in figure 3 is to be attributed to their very rapid 
switching-on. 

The charmed-quark contribution to i'W2 in muon scattering has been measured^^^^ and found 
to vary extremely rapidly with Q'^. It seems that this is a threshold effect, though we 
must emphasise that there is no good theoretical understanding of how one should take 
account of threshold effects^-*^^]. Such an understanding would require more knowledge than 
we have of the effects of confinement. Nevertheless, a successful phenomenology of the data 

IS 

available[llHl3]: 

the threshold effect is well described by supposing that the contribution 
to vW2 from charmed quarks is a function not of the Bjorken variable x but of 



Then 



^2 _ w2 



with Wq = + mp, and the best description of the data corresponds to choosing the 
parameter /ic such that Wq = niD + rn\^, the threshold value of W. (Strictly speaking, the 
threshold is different according to whether it is a quark or an antiquark that has absorbed the 
virtual photon, but the data are not sufficiently accurate for this to matter, and we simply 
take the lowest physical threshold.) Hence, by using the variable instead of x, we ensure 
that F2'^ goes to zero at threshold for each Q'^. According to the spectator-counting rulet^'^J 
F^^ should behave as {1 — as — 1, while Regge theory requires it to have power 
behaviour close to at small ^c- In addition, it must vanish linearly with when 

0, which we achieve as in (4a). This led us in reference 13 to the fit 



where we took e to be 0.086 and Cc — 0.045. For the reasons we have explained, we now 
adopt the value 0.0808 for e, but the difference is negligible. Since we made this fit, the data 
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have been scaled downwards, because of a revised experimental value of the branching ratio 
of the charmed quark decaying to a muon; we therefore shall use the value — .032. The 
fit, with the renormalised data, is shown in figure 4. 

Some discussion of the use of is called for. As we have indicated, it is suggested by the data. 
More naive theoretical considerations might lead to the "slow rescaling" choice jic = "t-c, the 
charmed quark mass, instead of the 2 GeV that is needed. To derive such slow rescaling one 
would need two assumptions: (i) that the mass scale associated with the fragmentation of the 
charmed quark after it has absorbed the virtual photon is rric rather than that of the lightest 
hadron to which it can fragment; and (ii) one can ignore the fact that the momentum k of 
the quark before it absorbs the photon is not on shell. The fact that needs to be so large 
indicates that neither of these assumptions is tenablet^l; not only is the minimum mass of the 
hadron to which the quark fragments important, so also is the minimum mass of the residual 
fragments of the proton left behind when the quark has been pulled out of it. The only way 
that such a mass can be supported is for the quark momentum k to go off shell and become 
negative: see (3). 

Given that the phenomenology of F^*^ points to the use of the variable defined in (5a), 
it is natural to assume that the strange antiquark distribution xs{x, Q^) should be handled 
similarly. The corresponding variable should be defined similarly, but with a scale jjis such 
that the threshold Wq — uik + wa- From (5b), this requires /x? — 1.7 GeV^. We shall work 
with 

xs{x,Q^) = C-^^y'\l-^sy (6) 

The threshold for the strange quark distribution xs{x,Q'^) is somewhat higher, but we shall 
not take account of this refinement, because we also have to ignore the possibility that the 
s and s quarks recombine^^^l and do not appear in the final state, with then a much smaller 
value for Wq and hence for Hs- We can only guess what to take for the parameter a^; we 
choose the value 1 GeV^, which is in between the corresponding value for light quarks (see 
(4)) and charmed quarks. Fortunately, its precise value will not be too important for us. Our 
first guess for the value of the constant C is that it is the same as for light quarks; (4) then 
gives C 0.22. We shall confirm from our later fit that this is the appropriate value in the 
case of the light quarks, but of course it is far from obvious that we should use it also for 
strange quarks. There is some confusion about the magnitude of the strange-quark content 
of the proton. A direct measurement from two- muon events in neutrino scattering t-*^^! finds 
that the strange quark distribution is half as large as the light antiquark distributions, but 
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the less direct methodl^^] of measuring 



xs{x, Q') ^ IF!^^ - 2,F^'' (7) 

gives a result that is significantly largerl-*^^], even allowing for that fact that it is risky to 
measure a small quantity as the difference between two large ones, as in (7). As a contribution 
to this debate, we point out that the threshold effects, being sensitive to the minimum mass 
that can be produced in the final state, will vary according to how one measures xs{x, Q"^). In 
there are two contributions. In the first, which is Cabbibo suppressed, the strange quark 
absorbs the weak current and becomes a light quark, so that the threshold is Wq = m^r + mA 
corresponding to Hs = 0.8 GeV. In the second, the strange quark absorbs the weak current 
and becomes a charm quark, so that the threshold is Wq = uid + rn^ corresponding to Hs = 
2.9 GeV. The latter process is also the source of the dimuon events in neutrino interactions, 
with the same value of (a much larger value than was assumed by the experimentalistst-'-^] 
in the analysis of their data!). In F^^ the strange quark survives and Wq = uik + tua, 
with then /ig = 1.3 GeV. In the latter case there is, however, the possibility that the s and s 
quarks recombine with consequently a smaller value for fXg- As an additional complication, 
the requirement that ^^^2 vanishes as — > applies only to photon-exchange processes and 
not to processes involving the non-conserved weak current. To some extent one can take this 
into account I-'^'^] through PCAC, but only for light-quark distributions. 

Having decided, albeit somewhat tentatively, how to handle the contributions to i'W2 from 
the heavy quarks (fortunately they are small), we now turn to the u and d quarks and their 
antiquarks. With the threshold Wq = uIt^ + nip, //=0.53 GeV, and we use this to define a 
variable ^ similar to (5a). When — > 0, ^ ~ A^^/Si', so that in order to have i'W2 vanishing 
linearly with we replace the small-a; behaviour (4) with 

i^W2 ~ yir°-''°V(Q') + Be'^'^'m') (8a) 

where 
and 

Aa-i(//2)-0-0808 = 0.604 56" i(//2)0-4525 ^ (g^) 

This is the behaviour we want when a; or ^ is small; when ^ is close to 1 we impose instead 
the spectator-counting rulesl-*^^] to determine the powers of (1 — ^). However, we find that if 
we simply multiply the terms in (4a) by such powers, they spoil the good fit of figure 3 to the 
small-x data. This is because, while at the smallest x-value 0.008 of the data in the figure. 
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(1 — xY is close to 1, at the largest x-value of 0.07 it is already as small as 0.6, and the effect 
of using instead (1 — is even more marked. In order to overcome this problem, we shall 
use the simple sum of powers for values of ^ less than some fixed value ^q, which we leave as 
a free parameter. For ^ > we match to each term in (8a) expressions of the form 

const 

where the power A is fixed by requiring such a form to fit smoothly on to the simple power of 
^ at ^ = ^0, that is the two forms and their first derivatives are equal there. We shall use the 
same value of throughout, and find that the best fit requires it to be quite small, about 
0.07. 

Consider first the valence distribution xuv{x,Q'^). At small ^ it should have the power 
behaviour 0.4525, corresponding to oj/p exchange. As ^ — > 1 its behaviour should be (1 — ^)^, 
according to the spectator-counting rule[14]. Hence we take 

We fix (3u and in terms of by requiring these two forms to join smoothly at ^ = ^o, 
and then determine by imposing the number sum rule 

' f = 2 (10a) 
s 

The valence distribution xdv{x,Q'^) should have the same power behaviour 0.4525. The 
spectator-counting rule would make it also behave as(l — ^)^as^— >1; however, it is well 
known that this conflicts with the measured ratiot-*^^] uW^' /fW^^^ which indicates that in xdy 
the coefficient of (1 — 0^ very small. So we take 

xdv{x,Q^)^D{^)tl;{Q^) 

with f3d and Ad again being fixed in terms of B^ by joining the two forms smoothly and then 
Bd determined from 



For the nonvalence distributions, for each hght quark and antiquark we include a term that 
behaves as ^-o osos for smaU ^: 

The constants 7 and As are fixed in terms of C by joining these two forms smoothly at ^ = ^o, 
with the same ^0 as for the valence terms. There is no number sum rule for the nonvalence 
distributions, so C is a free parameter. In iyW2, (Ha) is multiplied by 10/9, that is the 
constant A in (8a) is lOC/9. 

When the two-component parton model was first formulated^^^ , the component that later 
came to be called nonvalencel^-*^] had no contribution from p, a;, /, 02 exchange. This is be- 
cause, at that time, the idea of exchange degeneracy was taken very seriously, not just for 
the Regge trajectories but also for their couplings. Since then, it has become clear that the 
degeneracy of the trajectories is satisfied very well but there is no basis for supposing that 
it extends also to the couplings!^]. So we include in each nonvalence distribution also a term 
behaving as ^O '^^^s small ^. We write the total such contribution to i'W2 as 



where again there is only one free parameter, B, after we have joined the two forms smoothly. 
Notice that we have not specified how this term divides among the quark fiavours: while by 
definition (because the term is nonvalence) u and u receive equal contributions, as do d and 
(i, it is likely [22] that u ^ d. Notice also that we have, somewhat arbitrarily, chosen the power 
(1 — ^)^; the spectator-counting rule gives no guidance here, but we need a power greater 
than 7 to prevent the nonvalence distribution from becoming negative for large values of ^, 
since while Bu is close to | and Bd is about half that, the best fit for B is at quite a small 
value, as we already saw in figure 3. Again somewhat arbitrarily, we have used the same 
factor ip{Q'^) to make all the ^0.4525 ^gj-jns vanish linearly with Q"^ at small Q'^. 

The terms we have included so far have the the property that, at each x, vW^ increases as 
increases, while the data show the opposite when x is not small. At large Q'^ this is a 
consequence of the perturbative evolution, but we are constructing a fit at Q^-values where 
perturbation theory is not applicable. We therefore include an additional "higher-twist" term 

This resembles a term we introduced some years ago^^l, when we indentified it as a contri- 
bution from the virtual photon being absorbed by a diquark within the proton, though we 
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do not necessarily adhere to this interpretation now. By making it vanish quadratically as 
X — s> 0, we have ensured that it does not contribute to the real-photon cross-section, while 
the power (1 — 0^ gives a better fit than would (1 — ^). This term has two free parameters, D 
and Qo: in addition to the C and B we already have. Our best fit to the data for < 10 
GeV^ corresponds to 

C = 0.220 B = 0.279 = 0.071 D = 15.88 Qq = 550 MeV (13a) 

and is shown in figure 5. The per data point is just less than 0.5. Again we have quoted 
the values of the parameters to high accuracy because the errors, particularly in D and Qq, 
are strongly correlated. For example, 

C = 0.220 B = 0.274 = 0.075 D = 36.74 Qq = 338 MeV 

gives almost exactly the same x^- The values we obtain for C and B are rather stable, 
particularly that for C. The values of the subsidiary parameters corresponding to the choices 
(13a) are 

Bu = 1.456 Bd = 0.772 = 0.683 = 0.760 A = 1.144 A^ = 0.457 (14) 

We emphasise that our nonperturbative fit is supposed only to apply for Q^ < 10. For larger 
values of Q^ it lies above the data at moderate and large x, allowing room for the perturbative 
evolution to take over at Q^ = 10, or at some smaller value. See figure 6. We remark also 
that our choice to fit the data up to Q^ = 10 is somewhat arbitrary. If we replace this with 
Q^ = 5 the values of the parameters hardly change: 

C = 0.213 B = 0.312 ^0 = 0.069 D = 15.88 Qq = 554 MeV 

and the per data point is reduced to 0.3. 

We have extended our fit to i'W2, though the data here are rather more uncertain. Not only 
is there a lack of basic knowledge on exactly how to make deuterium correctionsE^^^, but also 
the NMC data are in the course of being changed slightlyt^^l. We use the NMC data, as 
published[19], for lyW^/uW^. Our fit to uW^ uses the same values for and C as for vW2, 
but allows for a different B, corresponding io u ^ and different "higher twist" . We choose 
to work with the same mass scale Qq in the latter, but allow its normalisation D to be chosen 
by the fit to the data. So we have 2 new free parameters. For we end up with 

C = 0.220 5 = 0.169 ^0 = 0.071 D - 4.94 Qo = 550 MeV (136) 
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which gives the fit shown in figure 7. Once again, the data do not determine the "higher 
twist" term at all precisely: we could have used a very different value of D with Qq changed 
correspondingly. But it is clear that we need different terms for and z/VF^, not only for 

the "higher twist" but also for the nonvalence ^-o o^os term: for its coefficient is -0.454, 
while for it is -0.337, so that indeed u ^ d. 

We have succeeded in obtaining an excellent description of i'W2 and i^Wj^/i^Wl' in the range 
< < 10 with very few parameters. The form of our fit is simple, and motivated by 
theoretical principles. Our analysis should be regarded as complementary to the standard 
onest^^] which are mainly driven by considerations of perturbative QCD. We do not use 
any neutrino-scattering data to determine the free parameters. As we have explained, until 
is rather larger than the range with which we are concerned, the neutrino data are not 
expected to be related at all precisely to the photon-exchange data, beacuse they approach 
the limit ^ Q differently. On the other hand, we fit the photon-exchange data down 
to small Q^, even — 0. The urgent question still remains of how properly to combine 
the perturbative and nonperturbative analyses, such as has been attempted for example by 
Badelek and Kwiecinski^^J. 

We have already remarked that the extrapolation of our fit down to very small x is finally 
not very different from that obtained from (4). In figure 8 we show our curves plotted 
against x. If the HERA experiments find results for VW2 significantly larger at small x than 
our extrapolations, we claim that this will be a clear signal that they have discovered new 
physics. Of course, the hope is that they will discover the Lipatov pomeron. In relation 
to this, the question arises whether the effect of the soft pomeron, which is included in our 
curves, should be subtracted off from the dataP], with the Lipatov pomeron being identified 
with anything that may remain, or whether instead the Lipatov pomeron simply replaces 
the soft pomeron in the small-x behaviour of vW2- This is all part of the general study of 
the interface between perturbative and nonperturbative QCD, which will be so important at 
HERA. 

One of us (PVL ) is pleased to acknowledge vigorous discussions with Jean-Rene Cudell, Steve 
Ellis and Dieter Haidt. 



10 



References 



1 A Donnachie and P V Landshoff, Nuclear Physics B244 (1984) 322 

2 H Abramowicz, E M Levin, A Levy, and U Maor, Physics Letters B269 (1991) 465; G A 
Schuler and T Sjostrand, preprint CERN-TH-6796/93 

3 J C CoUins and P V Landshoff, Physics Letters B276 (1992) 196 

4 J Forshaw, P Harriman and P Sutton, in preparation 

5 L Gribov, E M Levin and M G Ryskin, Physics Reports 100 (1983) 1 

6 P V Landshoff, Proc XXVII Rencontre de Moriond, (Editions Frontieres, 1992, ed Tran 
Thanh Van) 

7 A Donnachie and P V Landshoff, Physics Letters B296 (1992) 227 

8 ZEUS coUaboration, Physics Letters B293 (1992) 456; HI coUaboration, Physics Letters 
B299 (1993) 385 

9 P V Landshoff and J C Polkinghorne, Physics Reports 5C (1972) 1 

10 NMC coUaboration, Physics Letters B295 (1992) 159 

11 BFP collaboration. Physical Review Letters 45 (1980) 1465; EMC coUaboration, Nuclear 
Physics B213 (1983) 31 

12 P V Landshoff and D M Scott, Nuclear Physics B131 (1977) 172; V Barone, M Genovese, 
N N Nikolaev, E Predazzi and B Z Zakharov, Physics Letters B268 (1991) 279 and B304 
(1993) 176 

13 A Donnachie and P V Landshoff, Physics Letters B207 (1988) 319 

14 J Gunion, Physical Review DIO (1974) 242; R Blankenbecler and S J Brodsky, Physical 
Review DIO (1974) 2973 

15 A Donnachie and P V Landshoff, Nuclear Physics B112 (1976) 233 

16 CCFR collaboration. Physical Review Letters 70 (1993) 134 

17 CTEQ coUaboration, preprint FERMILAB-PUB-92-371 

18 A.D.Martin, W.J.Stirling and R.G.Roberts, Durham preprint DTP/93/22, RAL preprint 
RAL-93-027 (Proc Durham Workshop on HERA Physics, J Phys G, to be published) 

19 NMC collaboration. Nuclear Physics NP B371 (1992) 3 

20 P V Landshoff and J C Polkinghorne, Nuclear Physics B28 (1971) 225 

21 J Kuti and V F Weisskopf, Physical Review D4 (1971) 3418 

22 A D Martin, R G Roberts and W J Stirling, preprint RAL-92-078, Physics Letters B (in 

11 



press) 

23 A Donnachie and P V Landshoff, Physics Letters 95B (1980) 437 

24 P V Landshoff and J C Polkinghorne, Physical Review D18 (1978) 153; B Badelek 
and J Kwiecinski, Nuclear Physics B370 (1992) 278; M Strikman, XXVI International 
Conference on High Energy Physics (1992) 

25 B Badelek, talk at Durham Workshop, March 1993 

26 B Badelek and J Kwiecinski, Physics Letters B295 (1992) 263 



12 



Figure captions 

1 Data for the real-photon cross-section, with the fit from reference 7. 

2 The parton model 

3 NMC data at small x, with simple-power fit. Reading from top to bottom, the values of x 
are 

0.008, 0.0125, 0.0175, 0.025, 0.035, 0.05, 0.07 
For clarity of presentation, the curves and data have been scaled by a different factor at each 
value of x; the scale factors are 

4, 3.2, 2.5, 2, 1.5, 1.2, 1 

4 EMC data for xc(x, Q^) with fit from reference 13. The data have been renormalised to take 
account of the revised branching ratio for a charm quark decaying to a muon. Reading from 
top to bottom, the values of x are 0.00422, 0.0075, 0.0133, 0.0237, 0.0422, 0.075, 0.133, 0.237 . 
For clarity of presentation, the curves and data have been scaled by a different factor at each 
value of x; the scale factors are 128, 64, 32, 16, 8, 4, 2, 1 respectively. 

5 Data from NMC, SLAC and BCDMS with fit described in the text. For clarity of presentation, 
the curves and data have been scaled by a different factor at each value of x. Reading from 
top to bottom, the values of x, with the scale factors in brackets, are 

(a) 0.008 (5), 0.0125 (4), 0.0175 (3.2), 0.025 (2.5), 0.035 (2), 0.05 (1.5), 0.07 (1) 

(b) 0.09 (7), 0.1 (5), 0.11 (3.5), 0.14 (2.5), 0.18 (2), 0.225 (1.5), 0.275 (1) 

(c) 0.35 (32), 0.45 (16), 0.5 (8), 0.55 (4), 0.65 (2), 0.75 (1) 

6 The curves of figure 6b extrapolated to larger Q^. Reading from top to bottom, the values 
of a; are 0.09, 0.18, 0.275, 0.45 and 0.65, with scale factors 8, 4, 2, 1, and 1. 

7 NMC data for z/VF|*/z/VF2 , with fit. The data and curve correspond to values of that vary 
with x, from (Q^) =0.4 GeV^ at the smallest x to 10.8 at the largest. 

8 Fits to i'W2 extrapolated to very small x for three values of Q^. The curves are the total, 
together with its valence, nonvalence and "higher-twist" components 
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